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Qh' We prove that for compatible weakly nonlocal Hamiltonian and symplectic oper- 

• ators, hierarchies of infinitely many commuting local symmetries and conservation 

laws can be generated under some easily verified conditions no matter whether the 
I generating Nijenhuis operators are weakly nonlocal or not. We construct a recursion 

operator of the two dimensional periodic Volterra chain from its Lax representation 
and prove that it is a Nijenhuis operator. Furthermore we show this system is a 
I (generalised) bi-Hamiltonian system. Rather surprisingly, the product of its weakly 

Ch ' nonlocal Hamiltonian and symplectic operators gives rise to the square of the re- 

^ . cursion operator. 
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1 Introduction 



Integrable nonlinear evolution equations possess many hidden properties such as infinitely 
many symmetries and conservation laws. These symmetries can be generated by so-called 
recursion operators [HIS], which map a symmetry to a new symmetry. All known recursion 
operators including operands [3J for nonlinear integrable equations in the 2 + 1-dimension 
are Nijenhuis operators. The important property of such an operator is to construct an 
abelian Lie algebra. This property was independently studied by Fuchssteiner |1] and 
Magri [5], where they named the operator hereditary symmetry. For example, the famous 
Korteweg-de Vries (KdV) equation 

= Uxxx ~t~ QuUx 

possesses a recursion operator 

^ = Dl + Au + 2uxD-\ 

where stands for the left inverse of D^. Thus this recursion operator is only defined 
on ImDx- It is a Nijenhuis operator and generates the KdV hierarchy 

Ut^=^^M, j = 0,l,2,-- - . 

Any polynomial of 3? with constant coefficients such as 3?^ is also a recursion operator of 
the KdV. However, the operator 3?^ doesn't generate the whole KdV hierarchy starting 
form Use- 

The concept of Hamiltonian pairs was introduced by Magri [6]. He found that some sys- 
tems admitted two distinct but compatible Hamiltonian structures (Hamiltonian pairs) 
and named them twofold Hamiltonian system, nowadays known as bi-Hamiltonian sys- 
tems. The KdV equation is a bi-Hamiltonian system. It can be written 

ut = Dx 6{-^ + u^) = [Dl + AuDx + 2ux) 5{-), 

where 5 is variational derivative with respect to the dependent variable. These two dif- 
ferential operators and D'^ + AuDx + 'iux form a Hamiltonian pair. 

Interrelations between Hamiltonian pairs and Nijenhuis operators were discovered by 
Gel'fand & Dorfman [7] and Fuchssteiner & Fokas [SI E]. For example, the Nijenhuis 
recursion operator of the KdV equation can be obtained via the Hamiltonian pair, that 
is, 

^={Dl + AuDx + 2ux) D-^ . 

Such a decomposition of the operator 3? corresponds to the Lenard scheme used to con- 
struct the hierarchies of infinitely many symmetries and cosymmetries. The story why 
this concept was named after Lenard is told in |10j . 

In fact, the decomposition of 3? is not unique since it can also be represented as 

n=Dx {Dx + 2uD~' + 2D^\), 
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where the operator is Hamiltonian and the operator D^+^uD^ ^ + 2D^ is symplectic. 

The majority of 1 + 1- dimensional Hamiltonian integrable equations possess the same 
property as the KdV equation: their Nijenhuis recursion operators can be decomposed 
into the products of weakly nonlocal [TT] Hamiltonian and symplectic operators of order 
not less than —1 ( see p2] for a list of integrable systems). An exceptional recursion 
operator can be found in [13] although it can be represented as a ratio of two compatible 
weakly nonlocal Hamiltonian operators 



In applications, the Lenard scheme for both Hamiltonian and symplectic pairs [15] requires 
that one of the operators is invertible, which is not clearly defined for (pseudo-)differential 
operators in the sense that the inverses of many differential operators are no longer local. 
For compatible Hamiltonian and symplectic operators, there is no need to invert any op- 
erator in construction of Nijenhuis operators although the operators considered are likely 
to be nonlocal. The nonlocality has motivated Dorfman to introduce Dirac structures to 
the field of soliton theory [16]. Dirac structures unify and extend both Hamiltonian and 
symplectic structures. She showed that pairs of Dirac structures give rise to Nijenhuis 
relations, which is a generalization of the Nijenhuis operators associated with pairs of 
Hamiltonian structures, and further generalised the Lenard scheme [15] . 

In this paper, we restrict to weakly nonlocal differential operators. Without using Dirac 
structures, we prove that for compatible Hamiltonian and symplectic operators, hierar- 
chies of infinitely many commuting local symmetries and conservation laws can be gen- 
erated under some easily verified conditions. Their nonlocal terms suggest the starting 
points of the Lenard scheme. We treat the case when the generating Nijenhuis opera- 
tors are weakly nonlocal in Theorem [T] in section 13.11 and the case when the generating 
Nijenhuis operators may not be weakly nonlocal in Theorem [2] in section 13. 2[ 

We apply these results to study the algebraic structures of the following system 

Pi,t = (pLcccc + 2(^2,.. + 20i,,02,. + + 3e2<^^ - 3e-2*i-2^2 
p2,t = -201,.. - 02,.. - 20i,.02,. - 0i,. - Se^-^i + 3e-2'^i-2'^2 " ^ > 

We construct a recursion operator of system ([T]) from its Lax representation by apply- 
ing the idea in [17] for Lax pairs invariant under reduction groups [181 US |20l |2T] and 
prove it is a Nijenhuis operator in section 14. 2[ In section \^ we show this system is a bi- 
Hamiltonian system by constructing a Hamiltonian operator and a symplectic operator. 
Rather surprisingly, the product of its weakly nonlocal Hamiltonian and symplectic oper- 
ators gives rise to the square of the Nijenhuis recursion operator. The Nijenhuis operator 
itself does not possess such a decomposition. This phenomenon is very rare. The only 
known example to me in the literature is the system of gas dynamics which is of 
hydrodynamical type and where the operators are local. 

System ([T]) corresponds to the two dimensional periodic Volterra chain 

< n n _n 0„+7V = 0n, 2^ 0n = 0. (2) 

I ^n+l — + (Pn+l,x + (Pn,x — U, 

^ n=l 

with period = 3 [HI |19] . It has also appeared in the classification of integrable systems 
of nonlinear Schrodinger type |23j . 
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The two dimensional Volterra system ([2]) can be viewed as a discretisation of the Kadomtsev- 
Petviashvih equation. Indeed, in the hmit N oo, 

r = h^t, ^ = nh + Aht, rj = h^x 

system ((2]) goes to 

2 

O 

For integrable equations in the 2 + 1-dimension, their recursion operators are no longer 
pseudo-differential operators [211 ESI E]- To study the family of discrete integrable system 
(12]) for fixed N, we wish to shed some light on this issue. The exact solutions of system ([2]) 
have much in common with 2 + 1-dimensional integrable equations, which have recently 
been investigated by Bury and Mikhailov |26j . 

2 Definitions of geometric operators 

In this section, we sketch the basic definitions of Hamiltonian, symplectic and Nijenhuis 
operators following [TJ [151 122] • We begin with the construction of a complex of variational 
calculus. 

2.1 Complex of variational calculus 

Let X, t be the independent variables and m be a A^-dimensional vector-valued dependent 
variable, where G N is finite. All smooth functions depending on u and x-derivatives of 
u up to some finite, but unspecified order form a differential ring A with total x-derivation 

D^ = '^ Uk+i ■ ^ — , where Uk = d^u. 

fc=0 ^ 

Here ■ denotes the inner product of vectors. The highest order of x-derivative we refer to 
the order of a given function. For any element g & A, v^e define an equivalence class (or 
a functional) Jg hj saying that g and h are equivalent if and only if g — h & ImD^,. The 
space of functional, denoted by A', does not inherit the ring structure from A. 

Any derivation d on the ring A can be written as Yl'h=o^>' ' af^' 'where is an A^- 
dimensional vector field with entries from A. We denote the space of such vector fields as 
hk by A'^ . The derivation d is uniquely defined by its action on the dependent variable 
u and its x-derivatives. The derivation commuting with can be recovered from its 
action on the dependent variable, that is, Hq since we have hk = -D^/iq- This is known as 
an evolutionary vector field. Let f) denote the space of all such ho. For any P G f), there is 
a unique derivation dp = J2T=o ' af^- The natural commutator of derivations leads 
to the Lie bracket on I): 

[P, Q] = Dq[P]-Dp[Q], P,Qev, (3) 
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where Dq = Y^^o'^ ' -^x is the Frechet derivative of Q. 

We define the action of any element P G f) on j g E A! a,s follows: 

p / ^7 = / ap(<?) = / Er=o ^ ■ & = / [P] ■ (4) 

Direct computation shows that such an action is a representation of the Lie algebra f). 
Having a representation space of Lie algebra f), we can build an associated Lie algebra 
complex. This complex is called the complex of variational calculus. Let us give the first 
few steps since we do not need the general theory. 

We denote the space of functional n- forms by starting with = A! . Now we consider 
the space Q}. For any vertical 1-form on the ring A, i.e., uj = YlT=o^^ ' ^"^k, where 
h'' eA^, there is a natural non-degenerate pairing with the derivations dp: 

<u;,P>=J Er=o ■ D^^P = J {ET=oi-D.n') -P. (5) 

This can be viewed as the pairing of 1-forms of the form ^ du with P G P). Thus any element 
of is completely defined by ^ G A^ . For a given u, we have ^ = YlT=o(~-^^)''^'' ■ 

The pairing between Lie algebra i) and 1-forms allows us to give the definition of 
(formal) adjoint operators to linear (pseudo)-difTerential operators |28j . 

Definition 1. Given a linear operator S : i) ^ , we call the operator S* : [) ^ the 
adjoint operator of S if < SPi, P2 >=< S*P2, Pi >, where Pj G () for i = I, 2. 

Similarly, we can define the adjoint operator for an operator mapping from to i), from 
[) to P) or from ^l^ to Vl^. 

The variational derivative associates with each functional jg E A' its Euler-Lagrange 
expression S{fg) G defined so that 

< 6ifg), P >= {dfgm = Pfg =< Er=o(-^^)'&, ^ > > (6) 

where d : Q"' ^ Q^~^^ is a coboundary operator. Due to the non- degeneracy of the pairing 
(E]), we have 6{j g) = Yl'k'=oi~P^x)^'§^ G fi^. In the literature one often uses E referring 
to the Euler operator instead of S. For any ^ G fi^, by direct calculation we obtain 
d$, = — D^. We say that the 1-form ^ is closed if d^ = 0. 

Finally, we give the formulas of Lie derivatives along any X G f) using Frechet derivatives, 
cf. [15] for the details. 

Definition 2. Let Lk denote the Lie derivative along K We have 

LKjg = jDg[K] for JgeA'; 

Lxh = [K, h] for h E [); 

Lj^ = D^[K]+D*j,iO for ^ e fi^ 

LiSi = D^[K] - Dk^ + ^Dk for 3? : ^ [); 

LkH = Dn[K] - DkH - HD]^ for H : ^ i); 

LkI = Dx[K] + P)^X + IDk for I : [) n\ 
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In this complex we can identify all the important concepts in the study of integrable 
systems such as symmetries, cosymmetries, conservation laws and recursion operators. 
They are all characterised by the vanishing of the Lie derivatives with respect to a given 
evolution equation. This will be discussed further in section 14. 1[ 

2.2 Symplectic, Hamiltonian and Nijenhuis operators 

Definition 3. A linear operator S : i) —* (or ^ i)) is anti- symmetric if S = —S*. 

Given an anti-symmetric operator I : I) ^ Q^, there is an anti-symmetric 2-form associ- 
ated with it. Namely, 

uj{P,Q)=<IiP),Q>=- <IiQ),P>=-ujiQ,P), P,Qei). (7) 

Here the functional 2-form u has the canonical form [28] 

uj = ^jduAldu. (8) 

Definition 4. An operator X : i) is called symplectic if and only if the anti- 

symmetric 2-form ^ is closed, i.e., duj = 0. 

The symplecticity condition da; = can be presented in several explicit and equivalent 
versions. The details can be found in Theorem 6.1 in [T3] . 

For an anti-symmetric operator Ti. : —>■ i), we can define a Poisson bracket of two 
functionals 

{//, Jg}=<5if),n6ig)>. (9) 

Definition 5. The operator Ti is Hamiltonian if the Poisson bracket defined by ^ is 
anti- symmetric and satisfies the Jacobi identity 

Ih} + {{l9jh}, Jf} + {{JhJf}, Jg} = 0. 

For the Jacobi identity, there are several equivalent formulas given in [15] (see Theorem 
5.1). In [2H] (see p. 443), it was formulated as the vanishing of the functional tri- vector: 
J A D-hIHO] a 9 = 0, which makes it feasible to check. 

The Jacobi identity is a quadratic relation of the operator Ti. In general, the linear 
combination of two Hamiltonian operators is no longer Hamiltonian. If it is, we say 
that two such Hamiltonian operators form a Hamiltonian pair. Hamiltonian pairs play an 
important role in the theory of integrability. They naturally generate Nijenhuis operators. 

Definition 6. A linear operator 3? : f) P) zs called a Nijenhuis operator if it satisfies 
[^P, 3fJQ] - 3fJ[3fJP, Q] - 3fJ[P, 3fJQ] + ^^[P, Q] = 0, P,Q e f). (10) 
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Using the definition of Lie bracket ([3]), formula (fTOj) is equivalent to 

Lsfip^ = ^Lp^. (11) 

An equivalent formulation is: Dsr[5RF]((5) — ^Dsfi[P]{Q) is symmetric with respect to P 
and Q, cf. [9]. 

The properties of Nijenhuis operators [15] provide us with the explanation how the in- 
finitely many commuting symmetries and conservation laws of integrable equations arise. 



3 Lenard Scheme of Integrability 

The Lenard scheme was first used to generate the KdV hierarchy [10] . After the discovery 
of the interrelations between Hamiltonian pairs and Nijenhuis operators [71 [H] it was 
applied to bi-Hamiltonian systems. In 1987, Dorfman ([16j) introduced the concept of 
Dirac structures into the field of soliton theory in order to deal with nonlocal terms in 
the operators. 

In this section, we consider weakly nonlocal [llj Hamiltonian and symplectic operators, 
i.e., pseudo-differential operators with only a finite number of nonlocal terms of the form 
P C?) D^^Q, where P and Q are in the Lie algebra f) for Hamiltonian operators and in the 
space of 1-forms for symplectic operators. We prove that for compatible weakly non- 
local Hamiltonian and symplectic operators, hierarchies of commuting local symmetries 
and cosymmetries can be generated under some easily verified conditions without using 
Dirac structures. This is independent of the generating Nijenhuis operators being weakly 
nonlocal or not. Their nonlocal terms suggest the starting points of the Lenard scheme. 

Definition 7. A Hamiltonian operator Ti. : ^ i) and a symplectic operator I : f) — >■ fi^ 
are compatible if^ = 'HI is a Nijenhuis operator. 

We assume without loss of generality the nonlocal terms of a weakly nonlocal Hamiltonian 
operator 7i : fi^ ^ [) are of the form [2U] 

T.7=ie]Pj®D~^Pj, where e^ G {-1,0,1} and P,- G f) (12) 

and those of a symplectic operator 1 : \) ^ Vt^ are of the form 

ELi hlk ® D-^-fk, where h G {-1,0, 1}, 7^ G and D^, = (13) 

with the convention that if ej = or = 0, we take Pj = or 7^ = 0. Here (S> denotes 
the matrix product of two vectors (A^ x 1 column matrices), producing a N x N matrix. 

For a given weakly nonlocal operator S, its highest power of is the order of an operator. 
We say that the operator S is degenerate if there exists a non-zero weakly nonlocal 
operator T such that ST = 0. Otherwise, we say that the operator S is non-degenerate. 

Notice that the pairing between Pj G f) and 7^ appears in the computation of TCT. It 
is important to determine whether the pairing is zero or not. The pairings being zero 
implies that operator TiX is again weakly nonlocal. 
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Lemma 1. For any Q G f) and ^ G if D^ = then Lq^ = S<C,, Q>. 

Proof. Indeed, for any P G f) we have 

<6<C, Q>,P>=< D^[P],Q> + <^,Dq[P] >=<DIQ + D*q^, P> 
=<D^Q + D*Q^, P>=<Lq^,P>. 

Since the pairing is non-degenerate, we obtain Lq^ = 6<C,, Q>. ■ 

Proposition 1. Let the nonlocal terms of operators, Ti and I be of the form as [W\) and 
( flgj) . Assume that Pj, j = 1, ■ ■ ■ ,m and ^k, k = 1, ■ ■ ■ ,n are linear independent over 
C, respectively. If Lp.I = Lp^Ti = 0, then there exists anti-symmetric constant m x m 
matrix A^-'^ and n x n matrix B^^'' such that 

m n 

eiLp^P, = Y,PkA^i and iiLp^-fi = J^lkPl^i ■ 

k=l k=l 

Proof. The assumption Lp.I = Lp^Ti = imphes that 

ZT=i i^^Lp,P^ ® D~'P, + P,(g) D-h,Lp^P,) = 0; 
X;r=i [^iLp^li ® D-^-^i + 7i D-^eiLp^'ji) = 0. 

Applying Theorem [HI we obtain the resuhs as stated. ■ 

If for all J = 1, ■ ■ ■ ,m the matrices B^^^ are zero, the operator T-CT is weakly nonlocal 
since 6 < jijPj >= Lp.-ji = according to Lemma [H A lot of work has been done for 
this case, e.g. [301 Elj- We give the Lenard scheme including all the starting points in 
section ISTTl If there exists B^^^ ^ 0, the operator HX is no longer weakly nonlocal. The 
locality in this case has not been answered so far. In section 13.21 we tackle this problem 
and work out concrete examples. 

3.1 Case I: operator HX is weakly nonlocal 

Theorem 1. Let Ti and 2 be compatible Hamiltonian and symplectic operators, as defined 
above. Assume that Lp.Pi = Lp.^k = Lp.I = Lp.Ti = L-^^j.7s = and ITC^jk is closed, 
where j,l = 1, - ■ ■ ,m and k, s = 1, ■ ■ ■ ,n. Then for all i,ii = 0,1,2, ■■■ , 

1. = {XTCY'yk G are closed 1-forms and h\ = TtCk ^ commute; 

2. = {TiXyPj G f) commute and Q = Xp* G are closed 1-forms; 

3. vector fields h\ and pJ commute for all j = 1, ■ ■ ■ ,m, k = 1, - ■ ■ ,n. 

If there exist f I and gj such that = 6fl and Q = ^9), then all andp'j are Hamiltonian 
vector fields and their Hamiltonian are in involution. 
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Before we proceed with the proof, we first give a few lemmas. Part of Lemma [2] and 
Lemma [H] have been formulated and proved in abstract manner in [15] (see Proposition 
2.4 and 2.8). Here we give a straightforward proof. 

Lemma 2. Let I be a symplectic operator. For P G [) such that IP G , then IP is 
closed if and only if Lpl = 0. 

Proof. We know Djp = Djp if and only if for any Q,H & [), 

< Djp[Q],H>=< Dip{H),Q>=< Djp[H],Q>, 

that is, 

=< Djp[Q],H >-< Djp[H],Q > 

= < Dj[Q]{P),H > + < IDp[Q],H >-< Dj[H]{P),Q > - < IDp[H],Q > 
= -< Dj[P]{H),Q> - < D*pI[H],Q> - <IDp[HiQ>, 

where we used the fact that X is a symplectic operator. This leads to 

Dj[P] + D*pl + IDp = 0. 

From Definition [21 it follows that Lpl = 0. ■ 

Lemma 3. Let Ti he a Hamiltonian operator. IfHS, = P for some ^ E and = D^, 
then LpTi = 0. 

Proof. We know that LpH = D-^lP] — DpH — HDp. For any p,q E ^l^, we compute 
< {Lpn){p),q>=< Dnm]{p),q> - < Dndnp],q> - < nD^T,^{p),q> 

= -< Dn[nqm,p >-< nD^[np],q> + < Dndnq],p >= 0, 

where we used = and Ti being a Hamiltonian operator. This leads to LpTi. = 0. ■ 

Proof of Theorem [H The assumption that Lp.^k = and the closedness of ^ leads to 
S <7fc, Pj>= 0. Thus the operator 3? = HI is weakly nonlocal with nonlocal terms as 
follows: 

Besides, 3? is Nijenhuis since Ti and I are compatible. 

We first check the conditions of statement [2] in Lemma [7] in the Appendix. For the first 
statement in the theorem, we only need to check Lp^3fJ = L-^-yJ^ = 0. We know 

= Lp^ (HI) = Lp^ {n)I + HLp^ (X) = 0. 

Now we show L-yi^JR, = 0. Since Lp.^^ = 0, so 7^7^ is local. Together with D^^ = D*^, 
we have Ly^^^Ti = by Lemma El The assumption that ITi^k, which is local due to 
L'H^i^'ji = 0, is closed leads to L-n-yf.I = from Lemma El So L-n-^^Jk = Ln-^f^ilHI) = 0. 
Thus we prove that 1-forms = (IHY^k are local and closed. 
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For the second statement in the theorem, we need to show Lp^ilPi) = L-}i^^{TPj) = 
and TTCIPj is closed. 

It is clear that Lp^iTPi) = Lp.(I){Pi) +ILp.{Pi) = by the assumptions. We also know 

Lp^{n-/k) = Lp^Hi^k) + n{Lp^^k) = 0. (14) 

So Lq-i^^iTPj) = L-^^^(X)(Pj) + ILfi-y^iPj) = 0. From Lemma O IHTPj being closed 
is equivalent to Lt-ixp.I = L^p-I = 0. It follows from Lemma [2] and Lp.I = that 
XPj is closed. We also know ^Pj is local since we have proved Lp.{IPi) = 0. Thus 
Lsnp^Ti. = according to Lemma [31 Since 3? is Nijenhuis, we have Lsf^pM = ^Lp.^ = 0, 
that is, = {L^p.7i)T = TiL^p.T = TiL^p.T implying L^p.T = 0. We can now draw the 
conclusion that 1-forms Q = ^P] are local and closed. 

Next the conditions of statement [1] in Lemma [7] in Appendix are satisfied when we take 
into account what we have proved. Thus we can conclude for fixed k and j, vector fields 
hi = H^k commute and = {HJyPj commute. 

The third statement in the theorem follows from the fact 3? is Nijenhuis and 

Lp^Pi = Ln^,{n%) = Lp^in-fk) = 
for all J, I = 1, ■ ■ ■ ,m and k, s = 1, ■ ■ ■ ,n. 

Finally, we prove {fl, g^} = if and exist. The other cases {/^, /^^} = and 
{dji^ df} = can be proved in the same way. We have 

{fl gf} =< a^-HQ >=< {in)%, {nir^^p, >=< 7., ^'^^'^'p, > 

and 

6 < 7fc, 3fJ^+ip, >= L^s+rp^jk = 5 < in^k, 3fJ^P, >= L^sp^ilHjk) 
= {L^sp^{IH))-ik+1'HL^sp^{'-fk) =XHL^sp^{-fk). 

By induction, we obtain L^s+ip.-^f. = 6 < ■jk, ^'^'^^Pj >= implying |/^, g^^^ = 0. By 
now, we complete the proof. ■ 

This theorem gives rise to the Lenard scheme as shown in Figure [H for fixed k and j. 

In fact, such a scheme has been implicitly used for some integrable equations including 
the new systems in [3^. Since we have, for P, Q E I) or P, Q E , 

P ® D-'Q + Q0 D-'P = i(P + g) ® D-\P + g) - 1(P - Q) ® D-\P - Q), 

we can easily adapt Theorem [T] in terms of P and Q instead of P + Q and P — Q since 
all the operations involved in the theorem such as Lie derivatives and Poisson bracket are 
linear. 

Example 1. The Sawada-Kotera equation 

Ut = Ur^x + ^UU^x + ^U^U2x + Sw^Ma; (15) 
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Figure 1: Lenard scheme for compatible Hamiltonian and symplectic 
operators when the generating Nijenhuis operator is weakly nonlocal 

has a Hamiltonian operator 

n = Dl + 2uD^ + 2D^u 
and a compatible symplectic operator f3^ 

X = Dl + uD.^ + D.^u + {u^^ + y)^;^ + D-^{u.^x + —)■ 

Starting from both 1 and u^x + ^■^ '^i '^2 Figure [21 we generate the hierarchies 
of symmetries and cosymmetries of the Sawada-Kotera equation. 



3.2 Case II: operator TLX is not weakly nonlocal 

We start with a known example: although both the Hamiltonian and the symplectic 
operator are weakly nonlocal, the operator TiT is not weakly nonlocal. 

Example 2. The two-component system 



Ut = Uxxx + QuUxV + 3u Vx 
Vt = Vxxx + 9uvVx + Sv'^Ux 

possesses a Hamiltonian operator 

( -uD-'u Dx + uD-^v \ _ f Dx \ _f u . 
^ ~ V Dx + vD-'u -vD-'v J ~ \ Dx J y-vl^-^'' 

and a symplectic operator 

3vD-^v Dx + 2uD-^v + 3vD-^u 



(16) 



I 



Dx + 2vD-^u + 3uD~:^v 3uD-^u 



Dx 
Dx 



+ 3 (l) {v u) + D-^ {u v)- (_"^) {u -v) . 



These compatible Hamiltonian and symplectic structures of system fTB) first appeared in 
l^y, although the given symplectic operator was incorrect as stated there. 
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Using the notations in (fT2|) and (fT3l) . we have 




Indeed, 

/O -2 0\ 

(^Pi7i> -Lp^72, Lp^-fs) = (71, 72, 73) 2 7^ 

\0 0/ 

and thus we can not present TiX as a weakly nonlocal operator. 

Definition 8. We say Qi = Q2, where Qi,Q2 G f), with respect to {(3i,--- ,/9„}, where 
Pi G if Lq^^Q^pi = for i = 1, - ■ ■ ,n. 

Lemma 4. Let 7^, i = 1, 2, 3 and Pi be defined as [Tl\ ) in ExamplelM For any Q G f) and 
anti- symmetric 3x3 constant matrix A, if iiLq^i = Yl'k=i Ik^ki, then Q = Span^ < Pi > 
with respect to {71,72,73}- 

ti 

Proof. From Lq7i = A2172+A3173, we obtain Q = ^ {A21U + A-ii^/3v, A^i^/Su - A21V) 

with respect to 71. For such Q, it follows A32 = from —LqJ2 = — ^2i7i + ^3273 and 

tr 

further ^31 = from Lq73 = —743171. Thus we have Q = ^ (m, —v). ■ 

This lemma is inspired by Theorem [8] in Appendix. The idea is to identify Q G f) such 
that ikLq'jk = Yl^=i li^ik wheu the anti-symmetric n x n constant matrix A^ Q for given 
a set of linear independent 7^, i = 1, ■ ■ ■ ,n over C. Such an evolutionary vector field is 
rather strict. Therefore, we introduce the following definition. 

Definition 9. We say ^ E is proper for the operators Ti and I if for all 1 < I E N , 
vectors (7iX)'7i^ have no intersection with Span^ < Pi, - ■ ■ ,Pn>. 

Theorem 2. Let Hamiltonian and symplectic operators, Ti and X with nonlocal terms 
being (W^) and be compatible. Assume that 

1. Ti is non- degenerate; 

2. 7fc are linear independent over C for k = 1, ■ ■ ■ ,n; 

3. LpJ = Lp^n = for J = !,■■■ , m; 

4- For an anti- symmetric nxn constant matrix A and Q G P), if ikLq^k = Y17=i li^ik, 
then Q = Span^ < Pi, • ■ ■ , P„ > with respect to {71, ■ ■ ■ , 7„}. 

Then for a proper closed 1-form satisfying Lp.^^ = L^^^o'-yk = Lf^^oC,'^ = such that 
X7Y^° is closed, all = (X7Y)*^° G are closed 1-forms and W = 7^^* G f) commute for 
i = 0,1,2, ■■■ . If, moreover, = Sf\ then all W are Hamiltonian vector fields and their 
Hamiltonians are in involution. 
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Proof. Since Ti. and X are compatible, we have that 3? = TiX is Nijenhuis. From the 
definition of /i° and Lemma [2] and [H] this leads to LhoH = L^oX = and thus Lf^o^ = 0. 
Therefore, under the assumptions, if /i* G t) and & fl^, i.e., local, then the commute, 
Lfit^ = and the ^* are closed. Thus we only need to show that and ^* are local. 

Assume that h'-~^ and ^' are local and Lp^^'- = for / > 1. We show that /i' and are 
local and Lp^^''^^ = by induction. It follows from Lemma [1] that 6 < PjiE} >= 0. Thus 
we have = Ti^^ G f). Moreover, Li^iTi = from Lemma [31 Since 

= Lfii^ = (^Lfii'hCjl -\- 'HLj^iI = TtL/^iX, 

this leads to L/^iT = due to the non-degeneracy of Ti. This implies that 

ELi hiLf^nk ® D-'^k + 7fc ® D~'Lj,nk) = 0. (18) 

It follows from Theorem [H] in Appendix that iiL^i^i = J22=i^kAki, where the Aki are 
constant and Aki = —Aik- From assumption HI if the matrix A ^ 0, then contains the 
vector in Span^. < Pi, ■ ■ ■ ,Pn >, which cannot be true since is proper with respect 
to operator Ti and X. Therefore, the matrix A must be zero and this implies Lj^ij^ = 0. 
By Lemma [H we have 6 < h\^k >= and thus ^'"""^ = X/i' G VL^. Using the Leibnitz 
rule for the Lie derivative and the assumption Lp.X = Lp H = 0, we obtain Lp^'"*"^ = 

Lp^iine) = o. 

Finally, we prove that {/*, /*^} = if C,^ = 6f^ in the same way as we did for Theorem [H 
We have {f\ f'} =< i\nC' >=< iTnYi^ ,70X71)'' >=< C°,^''+'h^ > and 

= (LsRn+i-i^o(XH))i^° +X7YLsRn+i-i^o((^°) =I7iLy^n+i-i,^o{^°). 

By induction, we obtain Ljjii+i/jO^° = 5 < 3ft*i+*/i° >= implying {f\ /*^} = and 
this completes the proof. ■ 

This theorem gives rise to the Lenard scheme as shown in Figure O 




Figure 2: Lenard scheme for compatible Hamiltonian and symplectic 
operators when the generating Nijenhuis operator is not weakly nonlocal 

Before we apply it to concrete examples, we make a few remarks. 
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Remark 1. The assumption^^was inspired by Lemma^ From the proof of the theorem, 
the purpose of Definition O and this assumption is to enable us to draw the conclusion 
^^h^^k = from identity [W\) . For concrete examples, it is possible we do not require 
such strong assumptions. If is polynomial and we restrict Q to be polynomial, we only 
need to check that U does not contain linear terms in dependent variables since only such 
vector fields preserve the order and the degree ofjk- 

Remark 2. In application to nonlinear evolution equations, the 7^ in Theorem\^are good 
candidates for . If one of them, say 7„ without loss of generality, is indeed the starting 
point, we only need to check the assumptions\^and^in TheoremlEfor jk, k = 1, - ■ ■ , n—1. 

Proposition 2. Starting from = the hierarchy of commuting local symmetries 

and conservation laws for system [T^) (cf. Definition [77^) can be generated using the 
Lenard scheme as shown in Figure IE ■ 

Proof. To prove the statement, we check the conditions in Theorem [2l The operator Ti is 
non-degenerate since the determinant of the coefficient matrix of is non-zero. Following 
Remark [21 we only consider 71 and 72, which are obviously linear independent over C 
Secondly, by direct calculation we have LpJ-L = Lp^l = 0. Finally, as we remarked in 
Remark [1] we only need to check for all / that the vectors h'' do not contain any terms 
linear in u and v, which is true for the given operators and (although assumption H] has 
been proved in Lemma H]). 

Now we check the conditions on Notice that is closed and Lp-^^^ = 0. Moreover, 

So all the conditions of Theorem [2] are satisfied and thus we proved the statement. ■ 

Example 3. Consider the vector modified KdV equation JSBj 

Vt = V,,, + ?><V,V >V, + ?><V,V^> V, (19) 

where the dimension of vector V is N . Let Jij and Sij be N x N matrices. We represent 
its Hamiltonian operator and sympletic operator JWI as follows: 

n = D,+ {JijV) ^ D-\j,^vy, {j,j),i = 5^6'^ -sis'; 

l<i<j<N 

I = D, + 2V® D-'V + J2 -r-lJiS^JV) ® D-\S,,Vy, {S.,)u = 6^6] + 5\5l 

l<i<j<N- 

Proposition 3. Starting from = V the hierarchy of commuting local symmetries and 
conservation laws for equation ( fTPI) can be generated using the Lenard scheme as shown 
in FigurelE 

We first check assumption [H] of Theorem [2] in the following lemma. 

Lemma 5. The Lie derivatives ofTi andl as defined in Example\^ along the vector JijV 
vanish, that is, Lj v'H = Lj^.y! = 0. 
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Proof. According to Definition [TOl we have 
Lj.jv'H = Dn[JijV] — JijH + HJij 

l<k<l<N 

LjtjV^ = Dx[JijV] — JijI + XJij 

l<k<l<N 

Using the fact that 

[Jkli Jij] = Jkj^l + Jik^i + Jjl^k + Jli^k'^ 

[Skh Jij] = SkjS} — SikSj + Sji5l — Sii6l; 
we can prove that Lj-.yTi = Lj.^v^ = 0. ■ 

Proof of Proposition [31 We check the conditions in Theorem [2] one by one. First the 
operator Ti. is non-degenerate since the determinant of the coefficient matrix of is non- 
zero. For assumption [2], following Remark [2], we only consider SijV, 1 < ^ < j <, which 
are linear independent over C Assumption [3] is proved in Lemma [5l Finally, according 
to Remark [1], instead of checking assumption HI we only need to check for all / that the 
vectors /i' do not contain any terms linear in u and v, which is true for the given operators 
and 

Now we check the conditions on Notice that is closed and Lj-.yC^ = 0. Moreover, 
ne = V^; ine = V,, + 2<V,V>V = 5j -^{<V,,V^>-<V,V>^). 

So following the proof of Theorem [2] we obtain the results in the statement. ■ 

The following lemma show how to check assumption H] for Example 13 although it is not 
necessary since the objects we considered are differential polynomials. 

Lemma 6. Let JijV and SijV he defined as ^ in Example [H For any Q G f) and 
anti- symmetric constant matrix A, if 

^Lq{S,,V) = Y.i<k<i<N ^kuASuiV) , (20) 

where Aknj = -AijM, Akuj = Akij and Akuj = Auji, then Q = Span^KJij, l< i < j < N> 
with respect to {Sij,l < i < j < N}. 

Proof. When i = j, we know ^LqISuV) = Y^i<k<i<n'^kHt{SkiV). From Definition [2], it 
follows that the A;-th component of Q is equivalent to J2j'=i ^kjkkV^^\ where A^jkk = ^jkkk 
and y*^-'-' is the j-th. component of the vector V . For such Q and i < j, we compute 

Comparing to the coefficient of SuV in ( l20l) . we obtain that Amj = Ajijj = —Ajjij. This 
leads to Q = Y.i<i<j<n ^ijuJij^ with respect to {Sij, I <i < j < N}. ■ 
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4 Construction of recursion operators 



In this section, we construct a recursion operator of system ([T]) from its Lax representation. 
In general, it is not easy to construct a recursion operator for a given integrable equation 
although the explicit formula is given, cf. Definition [TOl The difficulty lies in how to 
determine the starting terms of 3?, i.e., the order of the operator, and how to construct 
its nonlocal terms. Many papers are devoted to this subject, see EHl EO]- If the 
Lax representation of the equation is known, there is an amazingly simple approach to 
construct a recursion operator proposed in [17]. The idea in [T7] can be developed for 
the Lax pairs that are invariant under the reduction groups. The general setting up and 
results will be published later. Here we only treat system ([T]). 



4.1 Integrability of evolution equations 

Before we proceed, we first give the basic definitions for symmetries, cosymmetries and 
recursion operators, etc. for evolution equations [ISl [2H] in the context of the variational 
complex described in Section [2l Meanwhile we fix the notation. 

To each element X G f), we can associate an evolution equation of the form 

ut = K. (21) 

Strictly speaking, this association is not as innocent as it looks, since one associates to 
the evolution equation the derivation 

1 + Td-k ^ 



m ^ ' dut 

k=0 

As long as objects concerned are time-independent as in this paper, one does not see the 
difference, but in the time dependent case one really has to treat the equation and its 
symmetries as living in different spaces, cf. [27] for details. 

Definition 10. Given an evolution equation f21\) . when Lie derivatives of the following 
vanish along K E i) we call: jg E A' a conserved density; h E i) a symmetry; ^ G 
a cosymmetry; 3? : f) ^ [) a recursion operator; a Hamiltonian operator Ti. : ^ i) a 
Hamiltonian operator for the equation; a symplectic operator X : \) ^ Vt^ a symplectic 
operator for the equation. 



Here we use the standard definition of a recursion operator in the literature. We refer the 
reader to ^9J for a discussion of the problems with this definition when symmetries are 
time-dependent. 

From the above definitions, we can show that if J/ is a conserved density of the equation, 
then 5{jf) is its cosymmetry. Moreover, if 7i is a Hamiltonian operator and T is a 
symplectic operator of a given equation, then Til is a recursion operator. Operator Ti 
maps cosymmetries to symmetries while 1 maps symmetries to cosymmetries. 
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We say that the evolution equation (12TI) is a Hamiltonian system if for a (pseudo-differential) 
Hamiltonian operator ?i, there exists a functional jf, called the Hamiltonian, such that 
Ti S(Jf) is a symmetry of the equation. Additionally, if for a (pseudo-differential) sym- 
plectic operator X, which is compatible with 7i, there exists a functional / g such that 

Iut=IK = 6{Jg) , 

we say that the evolutionary equation is a (generalised) bi-Hamiltonian system. 
The Sawada-Kotera equation (ITSl) is a bi-Hamiltonian system since we have 
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and 



4.2 Construction of a recursion operator 

Consider a matrix operator of the form 

L(A) = + A^VA - AA^V, (22) 
where A is the spectral parameter and A is a 3 x 3 matrix satisfying 

A = 

















and 












H 






^0 





.3^ 





e"^! 





e-^^ 










= A-1 = 


A2 



Notation 1. From now on, we often write Vi = e^' , i = 1, .., 3, where X]^=i 0i = 0. 

The operator L{X) is invariant under the following two transformations 

s : L(A) ^ 5-^L(aA),S and r : L(A) ^ -L*{\), 

A 

where S is a diagonal 3x3 matrix given by Su = a* and a = e^'^*/^. These two transfor- 
mations satisfy 

= = id, rsr = 

and therefore generate the dihedral group ©3. The reduction groups of Lax pairs have 
been studied in [I9l EDI El] • 

Assume that Y2^j=i ~ 0- Consider the zero curvature equation 

L(A), = [L(A), A(A)], (23) 



17 



where A{X) = X-^hVA - XA'^Vh + 3X-'^{VAf - 3X'^{A-^Vy and b is a diagonal 3x3 
matrix with entries b,, = (pi+i^x — (pi+2,x under the convention 0j = 0j mods for i > 3. It 
gives us a 3-component system 

(pi,t = (t)i+i,x. - (l)i+2,xx + - <Pi+2,x)(l)i,x + Se^'^^+' - 3e^'^^+' , 2 = 1,2,3. (24) 

Substituting the constraint J2'j=i 'i^j — i^^o ^^M^ "^^ obtain system ([T]). 

Notice that system ( l24l) is homogeneous if we assign the weights of Vi = e'^* as 1 and the 
weights of (pi as zero. We can also consider system ([T]) to be homogeneous under the same 
weights since we derive it from homogeneous system (1241) . This homogeneity enables us 
directly apply the results in 



The operator A{X) in ( |23l) is also invariant under the transformations s and r. In the 
commutator of the operators L{X) and A{X), the coefficients of positive powers of A are 
transposes of the negative powers and thus give no extra information. The constant term, 
i.e., the coefficient of A*^, is trivially satisfied. This is true in general. 

For given L{X), we can build up a hierarchy of nonlinear systems by choosing different 
operators ^(A) starting with A~"(l^A)". It is easy to check (VAY = I when J2'j=i 'Pj = 0. 
This implies that system ([T]) has no symmetries of order 3n. 

The idea to construct a recursion operator directly from a Lax representation is to relate 
the different operators A{X) using ansatz A{X) = VA{X) + R and then to find the relation 
between two fiows corresponding to A{X) and ^(A). Here V commutes with L{X) and R 
is the reminder \17\. 

Since the operator L{X) given by ( l22l) is invariant, we require that the ansatz VA{X) + R 
is also invariant. We take "P = A'^ + A~^, which is a primitive automorphic function of the 
group D3 and R is of the form 



i? = ^(A~Jrfj(FA)^ -A^(A-Vydj) where dj = \ dj,2 . (25) 
j=i \0 dj, 3 / 

This leads to 

L(A), = [VA{X) + R, L(A)] = VL{X)t + [R, L(A)]. (26) 

Substituting the ansatz (l25i) into (l26l) and collecting the coefficient of negative powers of 
A, we obtain 



A-^: VtA-VAd3 + d3VA = 0; 
A-3 
A-2 

A-i 



d2 - D.,ds - V^A4(yA)2 = 0; 

-A-^Vt - Dx{d2{VAy) - VAdiiVA) + A~^Vd^ + diiVAf - d^A'^V = 0; ^ ^ 



VrA = -Dx{diVA) + A-^Vd2{yAf - d2{yAfA-^V. 
We introduce the notation 
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and define A*0 = 0^*)A, where i = 1,2. Using this notation, we have Vt = V(pt and 
D^{VAy = ^^Zq cl)^J\VAy. Now formula can be simphfied as follows 

(l>t - 4'^ + ^^3 = 0; (28) 

D,d3 + 4'^ -4 = 0; (29) 

-2A-VtV^'A - D,{d2) - d2(p, - d2(pi^^ - 4'^ + di = 0; (30) 

(f)r = -D^{di) - + A-^d2l^'A - d2V^^^\ (31) 

Under the assumption Y^j=i 4'j = 0, we can solve system (l28l) . The unique solution for 
traceless matrix ^3 is 

rf3 = -^(2<^* + #) = ^(#-0*). 
Since (is is traceless, system ( l29l) is consistent. Its general solution is 

6?2 = h2D,d3 + Dj^^)+C2l =l^^2 + C2I, 



3' ' 3 



where C2 is a constant. In order to solve system flHOj) . the trace of its both sides should be 
equal, that is. 



= Tr(-2A-Vt^'A-D,(rf2)-t^20x-rf20i'^) 

1 

6 



-D,ie'^^ + e"^^ + e'^^) - 3D,C2 + -A(0L + 02. + 



where Tr denotes the trace of a matrix. So we can take 
We now substitute ^2 into flHUl) . It follows 



-20f V^(^) - |(/>i5, - (D.C2)/ + + C24^^ - 4^^ + 4 = 0. 

Solving for di, we obtain 

where Ci is a constant. Again due to the consistence of system flHTl) . ci satisfies 

= -3D^ci -Tr(di0x--4^^^^^^^ +4V'^^^^). 
Using the solutions of d2 and di, we can write it as 



SD^^Ci = ^ ( ai, ^2, as ) ( 02, 



i,t 
t 
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) {Dl — (j)i^xDx + 6vf) + 3 {vf_^_i — vfj^^Dx ■ Now we have determined 



where = 

^2 and di. Substituting them into (IHTj) . we have 

'Pi+l,xxt ~ '^4'i+l,x4'i+l,xt + C2<Pi+l,x + 3Ci) 
+ 'Pi,x {—2(l)i+2,tvf^2 ~ \4'i+2,xxt + \4'i+2,x4>i+2,xt — C2<pi+2,x 



i+l ' 3'+'i+l,xxt 



+ {4>i+l,xtvf^2 ~ 4'i+2,xtvf+i 



'^(Pi+l,x4'i+l,xt + C2<Pi+l,x + 3Ci) 



W+2J 



Dx {^—2(f)i+2,tvf+2 ~ \4'i+2,xxt + \4'i+2,x4'i+2,xt + 2(j)i^i^tvfj^i + |0i+l,a;a;t — \4'i+l,x4'i+l,xt) 
/o.,2^_ _ ^(f)-^^(f)-^^^{^cl)i^i^x ~ 4>i+2,x) + {4>i+l,xt'^'i+2 ~ 4'i+2,xt'vf^i) 



''j,x'rj,xt 



-| Ei=l((</'j+2,x - 0j+l,x)(6w|0j,t + - 

+ {~2(j)i+2,tvfj^2 ~ \4'i+2,xxt + \4'i+2,x4'i+2,xt + 20j+i_tuf_,_i + |0i+l,a;a;t ~ \4'i+l,x4'i+l,xt) 
-\ (3(ff+i - 1^,^+2) + {.4>i+l,xx - 4'i+2,xx) + 0i,a;(</>i+l,x " 0i+2,x)) Z]i=l 'i^j.^^'/'j,* 

Z)j=l((0i+2,x - - <i)j,x<Pj,t) + 3(wJ+i - 'y|+2)0i,t " {(Pj+2,xx " 

+ 1 (3(ff+i - ^,^^.2) + - (t)i+2,xx) + - 0i+2,x)) -D^^ Ej=l(2^^j + |0J>x)0j,t 

-\(t)i^xD~^ Yl]=l {{.4>j+2,x - (j)j+l,x){'ov'^j + 4>j,xx) + {,4>j+2,xxx " 4'j+l,xxx) 
+(t)j,x{(pj+2,xx - <Pj+l,xx) - 6(t;J+i0j+i,^. - t;|+20j+2,x)) • 



We substitute 03 = — 0i — 02 into the above expression. This leads to a recursion operator 
3? of system ([T]) mapping the flow ( 02,i ) to the flow ( 01,,-, 02,t ) , that is, 



^2.r 



3? 



4>i,t 

4>2,t 



3^21 



^12 
3f?22 



Theorem 3. System (Qp possesses a recursion operator of order 3 with entries 



3f? 
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3? 
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:0i,02) = ii^,' + i(0M + 02,.)i^^ 

+ |(301,xx - 01,x02,x + 302,xx - (Plx - <Plx + 21^^3 + + H)Dx 

-\iy4>l,x4>2,xx + 01,xx02,x + 20i,^0i,a;a; + 2(p2,x(p2,xx + 40L + 70? a,02,a' + 01,x<; 

+ i(5t;2 - 4t;i - 7t;|)0i,, + i(8t;2 - _ I3vl)(l>2,x 
-g4'i,xD~^S2,^^ - g0i,fL>^^Si,<^^; 

01, 02) = §1)3 + |0i,.D2 ^ 1(30^^^^ _ 202 ^, - 201,^02,. - 20|, + 6vj + 2Avl + 24i;|)D, 



^Ix) 



+ |(01,x01,xx - 01,x02,xx - (t)2,x4>l,xx " 202,x02,xx " 0?,^ " 4>l,x4>2,x " 01,a;0i,a;J 

+ §(^2 + vl- 5t;2)0i,^ + 4(t;| - v|)02,x - |0i,x^^^S2,</,2 - 101,4^-^51,^2; 



and 3^21 ( 



-3?12(02, 0l); 3^22(01,02 



<^^L,X^X "^,<P2 9" 

-3^11 (02, 01 ). Here 



51, ?ii 

52, </.i 



-201,.. - 02... - 6e2<^i + 6e-2<^^-2<^^ 
-01,.. - 202,.. - 6e2<^^ + 6e-2'^i-2<^2 



-302,... - 40i,.0i,.. + 202,.02,.. - 20i,..02,. - 20i,.02,; 



-12e2<^i 



n,x 



I2e'^^^2,x - 12e- 



2</.i-2</.2, 



'2,. 



S2,<^2 = 30 



201 ,.01,.. 402,.02,.. + 20i,..02,. + 20i,.0: 



xY2,xx 



+12e2'^i0i,. + 12e2<^20i,. + 126" 



■201-202, 



'1,. 



(32) 
(33) 



and 0i,t , 02,t is the system itself. 
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In Theorem [3l (132|) and (!33|) are two co-symmetries of system ([T]) . They are variational 
derivatives of the following two conservation laws: 

H2 = S(f)i^x4>2,xx + 4>l,x4>2,x - h,x<Pl,x + - |02,x + 6t'f(/>3,x + 6t;|(/)i,a; + Qv^4>2,x- 

The nonlocal terms of the recursion operator are determined by the symmetries and 
co-symmetries of the corresponding orders. Such structures of recursion operators have 
been discussed in [301 EI] • It has been shown that such recursion operator gives rise to 
hierarchies of infinitely many commuting local symmetries if it is a Njienhuis operator. 

Theorem 4. The operator^ defined in Theorem\^is Nijenhuis. 



Proof. We need to check that the expression H := Dy^[^P](Q) — 'SiDyi[P]{Q) is sym- 
metric with respect to two-component vectors P and Q. We use subindex i to denote 
i-th component. The first component of H will be written as Hi. The calculation is 
straightforward, but rather complicated. Here we only pick out the constant terms in H, 
i.e., terms are independent of the dependent variables 4>2 and their x-derivatives. We 
denote these terms by H^. 

For the recursion operator 3?, its constant and linear terms are 

,_!_( Dl 2Dl 
- 3 I -2Dl -Dl 



3\ -<j)2,xDl - D^(f)2,xD^ -{(t>i,x + (f>2,x)Dl-D^{4>i^^ + (t)2,x)D.^ 

So = Dyii[^°P]{Q) - 3ft°D5Ri[P](g). Due to the relations among the entries of the 
operator 3ft, it is easy to see the second component of i.e., (i^°)2 is related to its first 
component as follows: 

{H%{Pi,P2;Qi,Q2) = {H^)i{P2,PuQ2,Qi). 

Thus we only require to check whether the first component is symmetric with respect to 
P and Q or not. Notice that 

9{H%{Pi,P2;Qi,Q2) 

= {P2,ix — Pl,ix)Ql,2x + Dx {{P2,ix — Pl,ix)Ql,x + {Pl,ix + '^P2,Ax)Q2,x) 
+ {Pl,ix + 2P2,4x)Q2,2a; — Dl{{Pl,x — P2,x)Ql,2x — {Pl,x + '^P2,x)Q2,2x) 
~P^x{{Pl,x + P2,x)Ql,x + Pl,xQ2,x ~ 2(Pl,x + P2,x)Q2,x ~ '^P2,xQl,x) 

= {P2,5x — Pl,5x)Ql,x + {Pl,5x + 2P2,5a;)Q2,x + {P2,x ~ Pl,x)Ql,5x + {Pl,x + '^P2,x)Q2,5x 
+3-D^ {{P2,2x — Pl,2x)Ql,2x + {Pl,2x + 2P2,2a;)Q2,2x) 
—3 {{P2,3x — Pl,3x)Ql,3x + {Pl,3x + 2^2,30; )Q2,3a:) 

^P'x iiPl,x + P2,x)Ql,x + Pl,xQ2,x — 2(Pl,x + P2,x)Q2,x ~ 2P2,xQl,x) • 

This is symmetric with respect to P and Q and thus we proved the statement. ■ 

This Nijenhuis operator 3ft defined in Theorem [3] has two seeds: the trivial symmetry 
Ux = { 4>i,x, 4'2,x ) and system ([T]). We can generate the local symmetries of order 
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symmetries cosymmetries 




Figure 3: Interrelations between 3? and 7Y,X of system ([T]) 

3A; + 1 and 3k + 2 for A; > of system ([T]) by recursively applying the operator 3? on these 
two seeds. In particular, ^{ux) gives us a local symmetry of order 4 with first component 

-U'I'lx + (Pl,x(f>2,x + (/>lx)i<Pl,2x + 202,2x) - ^0t, ' l4>lx^lx ' |0?,x'/'2,.. ' 
_e-2*2 _ e-4<^l-4</.2 + g20i+202 ^ g4<^2 ^ |(60L + 12(^i,,(^2,x + 201,2. + 402,2x)e2<^i 

-i(202 _ 40i,,.02,. - 130t - 201,2. - lO02,2x)e2'^^ 
-|(702_, + 2201,^02,. + 130|, - 801,2. - lO02,2.)e-2'^^-2<^^ 

and the second component G(0i,02)</.2 = "^(02, 0i)</.i- 

The adjoint operator of 3?* gives rise to the cosymmetries of order 3n + 2 and 3n + 3. In 
figure [3], we list out the orders of symmetries in the left row and orders of cosymmetries 
in the right row. We use a circle around a number k to denote that the system does not 
possess the symmetries or cosymmetries of order k. 



5 Symplectic and Hamiltonian structures 

In this section, we show system ([T]) is generalised bi-Hamiltonian by presenting its Hamil- 
tonian and symplectic operators. Surprisingly, the product of these two operator does not 
lead to the recursion operator we constructed in section 14.21 but to its square. 

We know Hamiltonian operators map cosymmetries to symmetries while symplectic op- 
erators map symmetries to cosymmetries. For system ([T]), from Figure [3] we can draw the 
conclusion that the possible order of Hamiltonian operators can only be 3A; -|- 2 and of 
symplectic operators 3A; -|- 1. Here we consider positive orders, i.e., > 0. 

For system ([T]), there exists an anti-symmetric operator Ti such that 

V ^21 ^22 / V -51,02 / V 3^21 3fC22 / V ^2,x / 
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where 5R, s\ and sf are defined in Theorem [3] and 

7Yi2 = 3i^^ + 20i,^D^ - 2(f)2,xDx + (l)l,xx - <p2,xx - 4>\x - 4>\x - 4>l,x(p2,x 

+3ivf + vl + vl)-(j)i,xD-^(j)2,t-<pi,tD-^<p2,x 

and 7i2i(0i, = -^12(^2, 4>i)] ^22(^1, ^2) = -^ii(02, 0i)- 
Theorem 5. T/ie operator Ti defined above is Hamiltonian. 

Proof. We prove this by checking the condition of Theorem 7.8 in [28j. The associated 
bi-vector of Ti is by definition 

e = ^JeAne, where 9= {61,62)'^. 

We need to check the vanishing of the tri-vector: PiVfig^Q) = 0. Instead of writing out 
the full calculation, we pick out terms with highest degree in x-derivatives of dependant 
variables (pi and (/)2 in 3- form of 9i. The relevant terms in TiO are 

{{2(pl,x(p2,x + 0?,J^l)_i - (201,.;02,x + ^ 
-<p2,x {{2(l)l,x(p2,x + <PlxWl)_^ + i2<Pi^x<p2,x + 0i,J {(j)l,xOl)_i J ' 

where denotes D~^{p) and such terms in G are 

/ -<Pl,xei A ((2,^i,,,^2,. + <PlM -1 = 1 ^i^h- ^ {{2<Pl,x<p2,x + (Pl.)0i) _^ 

= -J {2cf),^,(p2,x + (pl^Wi A 
using integration by parts. Thus the terms we look at in Prv-^e(O) are 

/ -{2<f)i,xh,x + {(Pl,x0l)_, A ^1,, A {{2(j),,xCj)2,x + (Plx)Oi)_, 

+ J2{(j)2,x + <Pl,x)(pl,xx {{2(pl,x<p2,x + <PlxWl)_^ A A {(/)l,xOl)_^ 
+ J 201,x.02,x... {{2(pl,x^2,x + <Plx)0l)_i A A i(pl,xdl)-i 
= 0. 

Similarly we can prove that the tri-vector Prv^6i(0) vanishes, which implies that 7i is a 
Hamiltonian operator. ■ 

We now construct a symplectic operator of system ([T]). Its lowest positive order is 1. 
Proposition 4. No weakly nonlocal symplectic operator of order 1 exists for system (J\). 

Proof. Assume that any of the nonlocal terms of the symplectic operator is of the form 
^iD~^^2, where G Q^. Since system ([T]) is homogeneous in the variables (pt and Vi, 
is also homogeneous; its possible weights are 0,1. From the results in [30], the are 
cosymmetries of the system. One can check that the system has no cosymmetries of 
weight 1 and 0, that is, no conserved densities of weight and 1. ■ 
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We now look for next lowest order symplectic operator, which is 4. Indeed, there exists 
an operator X such that 

where 

9 = 3(/)?^3^ + 3(/)i,3^(/)2,3x + 30^3^ - ^<\)i^xx<\>2,xx{<\>\,xx + <\>2,xx) " (24^;^ + 15t;| + 2Avl)(f)l^^ 
+5(02, + 0i,,02,x + 0|J(0?,.... + <l)i,xx<p2,xx + 0|.J - {15vf + 24t;i + 24t;|)0i^,, 

-{15vj + 15vl + 33t;i)0l,xx02,xx + f (2v? - - 

+ f (2t>| - Vl - Vl)4>l^(t)2,xx - l<Plx<Plx{(l>l,xx + h,xx) + |(0tx'/'2,xx + 01,xx(/>2,x) 

+^ + 0U + U^lxh,x + Mix) + U^vl + 0|J4. + UH + 

+f + + fivl + ^l)<t>ix - m.^4>ix + {^vi - 'ivl + 'ivDMl^ 

+ (f -I - f + 'ivDMl^ + (223t;! - hvl - 5vM,x<l>lx 
+ {15vl + 60vf + 60vl + 66vjvl + 66vlv1 + S4vfv1)(l)l^ 
+ {15vf + 60vl + 60vl + 66vlv^ + 66vfvl + 84t;|t;|)0|, 

+{15vf + 15v| + 105v| + 48^2^;! + SAvjvl + SAv^vl)^i^a:(l>2,x - ^{vj + v^ + v^f 

and 

^11 = (201,0, + <p2,x)Dl + Dl{2(pi^^ + 02,x) + qD^ + D^q + |(si,<^,D-is2,0i + S2,,^iI^^^Si,<^J 
Ti2 = 3L»^ + 2(01,, - (f>2,x)Dl + p2Dl +p^D,+po + \{si,^,D-^ S2,<^, + S2,^,D-isi,^J. 

Here Sj,<^^, Sj,<^2 the components of cosymmetries and 

q — — 201,,,, — 02,a;a;a; — '2,(j)i^x<p2,xx — 4>2,x4>2,xx — ^4>l,x ~ 20f^,02,a; + |02,a; 

+(1302,0. - ^i,x)vl - (401,, + 202,^)t;| - (401,, + 1702,^)t;|; 

P2 = 30i,o;o; - 302,a;a; - 4>l,x " 4>l,xh,x " 02,a; + 15(?J? + t)l + vf); 

Pi = <Pl,xxx - h,xxx - 20l,a;02,a;a; " 402,x02,a;a; " |0"l,x- - 4>l,x4>2,x + 01,x02,x + t4>2,x 

+(1101,^ + 402,^)t;? - (401,, - A9(p2,x)vl - (1901,^ + 4102,^)^;!; 

PO = -(01,x + 2(j)2,x)h,xxx + 20f^,, + 01,a;a;02,a;a; - 2(0?_, + 01,a;02,a;)01,xx " {(p^x " 20i,a;)02. 

+(1001,0.0. + 202,:.:. - 20?,, - - (20i,c.a. - 2902,:.:. + 40?,, + 801,0.02,0. - 360i,,) 

-(1401,:.:. + 2502,:.:. - 60?,, - 3801,0.02,:. - 280|,)t;32 + 12(vl + + vl)^; 

And 22l(01,02) = -2^12(02, 0l); 2:22(01,02) = -2^11(02, 0l). 

Theorem 6. The operator I defined above is symplectic. 

Proof. We use the notation d0 = (d0i, d02)"'". The 2-form defined by the operator I is 

u; = i J d0 A Jd0 

= / ((201,, + 02,o.)d0i A d0i,3o. + g(0i, 02) d0i A d0i,, + |si,^^d0i A L>^^(s2,0id0i) 

+3 d0i A d02,4x + 2(01,0. - 02,x)d0l A d02,3o: + P2d0l A d02,2o: + Pld0i A d02,o. 

+Pod0i A d02 + |si,0id0iD-^(s2,</,2d02) + |s2,^id0ii:'-i(si,^2d02) 

-(202,0. + 01,x)d02 A d02,3o: - ?(02, 0l) d02 A d02,o. + |si,</,2d02 A i:*"^ (s2,</,2d02)) • 

We now compute da;. Instead of carrying out the whole computation, we only demonstrate 
the method by picking out 3- forms of 0i and its x-derivatives in da;. These are as follows: 

/ (2d0i,o; A d0i A d0i,3o; - 2d0i,3o; A d(f)i A d0i,o; - |d0i,o.o. A d0i A i:)~^(s2,</,id0i) 

-|si,^,d0i A L>-i((40i,, + 202,o;)d0i,,, A d0i + (40i,„ + 202,o.o.)d0i,x- A d0i)) 
= / (|s2,^i£'-^(d0i,,, A d0i) A d0i - |si,^jd0i A (40i,, + 202,o.)d0i,, A d0i) 
= / |s2,</,id0i,o. A d0i A d0i = 0. 
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By working out for other terms, we can show du = 0. Thus we prove T is symplectic. ■ 
Theorem 7. System ^ is a hi-Hamiltonian system. 

Proof. We only need to show that the symplectic operator X defined in Theorem [6] and 
Hamiltonian operator Ti defined in Theorem [5] are compatible. Operator Ti is of order 
2 and operator X is of order 4. This leads to that TiX is of order 6. By straightforward 
computation, one can verifies that TiX = 273?^, where 3? is defined in Theorem [31 From 
Theorem m we know 3? is Nijenhuis. So is the operator 3?^. Thus these two operators are 
compatible. Hence we proved the statement. ■ 

We can check that the conditions in Theorem [1] are satisfied for 7i and X. Therefore, 
starting from four starting points, two symmetries appeared in Ti. and two cosymmetries 
appeared in X, we can generate a hierarchy of commuting local symmetries, which are all 
Hamiltonian vector fields and their Hamiltonian are in involution. 



6 Discussion 



In this paper, we prove that for compatible weakly nonlocal Hamiltonian and symplectic 
operators, hierarchies of infinitely many commuting local symmetries and conservation 
laws can be generated under some easily verified conditions no matter whether the gener- 
ating Nijenhuis operators are weakly nonlocal or not. The problem how to generate local 
symmetries and conservation laws when the Nijenhuis operators are no longer weakly 
nonlocal has not been studied before. As in Example [2] and [3] where the objects are differ- 
ential polynomials, we believe that assumption H] in Theorem [2] can be relaxed in general. 
However, we are not able to simplify this assumption yet. 

We construct a recursion operator 3?, a Hamiltonian operator 7i and a symplectic operator 
X for system ([T]). We show that Til = 3?^. This leads to '^^^Ti being Hamiltonian and 
compatible to Ti for A; G N. An immediate question is: is 3ft7i Hamiltonian and compatible 
to HI We conjecture the answer is positive. However, the computation involved is rather 
big and we have not found an elegant way to prove it. 

In the Lax representation L{\) of system ([T]), cf. formula (1221) . A is a 3 x 3 matrix. The 
natural generalisation is A being n x n matrix. The construction in this paper works for 
any finite n. For a given n, the corresponding system possesses a recursion operator 3? of 
order n, which can be constructed in the same manner as in Section 14.21 The system is 
bi-Hamiltonian with the lowest positive order of Hamiltonian operator Ti. being n — 1 and 
that of a symplectic operator X being n + 1. These operators have the same properties as 
we discussed for n = 3, namely, HT does not give rise to 3?, but 3?^. 

If we treat arbitrary n by considering A as a shift operator, this leads to a 2+1- dimensional 
lattice-field integrable equation. Recently, Blaszak and Szum have constructed Hamilto- 
nian operators for such type of equations with a certain type of Lax operator [IQ]. It 
would be interesting to construct the bi-Hamiltonian structure for the 2 + 1-dimensional 
lattice- field equation and to see how this structure is related to the ones with finite periods. 
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Appendix 

Here we give the lemmas used in proof of Theorem [1] and Theorem O Lemma [7] is due to 
Sergyeyev [31]. To increase the readability, we include the statement here. 

Lemma 7. Let the nonlocal terms of a Nijenhuis operator 3? : f) ^ f) be of the form 



1. If for all j, k = 1, - ■ ■ ,1, we have Lq lSk = and both (3j and 3?*(/3j) are closed, then 
for any /iq G t) such that L^o'^ft = and Lhgf3j = 0, all hi = 3?*(/io) are local and 
commute, where z = 0, 1, 2, ■ ■ ■ . 

2. If Lq.^ = 0, then for any such that Lq^^q = and both and ^*{Co) o^re closed, 
all = 3?**((^o) fl^e local and closed, where i = 0,1,2, ■■■ . 

The following theorem was proved by A.V. Mikhailov, Applied Mathematics Department, 
University of Leeds and published here with his kind permission. 

Theorem 8. Let ffc,gfc G (either in f) and/or in Q^), vector- columns fi, ... , f „ be 

linearly independent over C and 



Y!j=i Qj ® D-^Pj^ where Pj E f) and (3j E n^. 



n 




(36) 



k=l 



then 



n 




i=l 



where 



0. 



The proof of the Theorem is based on two Lemmas. 



Lemma 8. // 



n 







k=l 



then for any p,q E 



n 




(37) 



k=l 
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Proof. Indeed, it follows from (l36l) that 



j=Q 



and thus for any q G Z>o 

n 

Y,ikDl{gl') + g,Dl{it) = 0, q = 0,l, 



{31 



k=l 



In order to demonstrate (IH71) we use induction. It follows from flHSl) that the statement is 
true for p = and any q G Z>o. Let us assume that ( 1371) is valid for any q and p < I — I 
and then show that it is true for p = I. For p = / — 1 we have 



k=l 



Applying we get 



k=l k=l 

The last sum vanishes due to the induction assumption. Thus (1371) is true for p = 
Let Ffc, Gfc denote infinite dimensional vector-columns 



/ ffc \ 



Dm 



\ 



/ gfe \ 

D.igk) 

msk) 



J 



and F, G denote matrices 

F=(Fi, ...,F„), 
Then flH7|) can be written in the form 



G — (Gi, . . . , G„). 



FG*^ + GF^^ = . 



(39) 



Lemma 9. Let vectors fi, . . . , f„ be linearly independent over C, then vectors Fi, . . . , F„ 
are linearly independent over A and thus 

rank (F) = n. 

Proof. Let us assume the opposite, i.e. rank (F) = m < n. Without a loss of generality 
we shall assume that the first m vectors Fi, . . . , F^, are linearly independent over A and 
thus the rest vectors F^, k = m + 1, . . .n can be expressed as 



^F^ask, ask ^ A, k = m + 1, . . . ,n. 



s=l 
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The latter is equivalent to 

m 

^^(ffc) = I]^^(fs)«sfe, fc = m + l,...,n, pGZ>o. (40) 
Differentiating (HOll we find 

m mm 
s=l s=l s=l 

and thus 

m 

s=l 

Since vectors Fi, . . . , are assumed to be linearly independent, we have D^^agk) = 
and thus ask £ C. It follows from (HOl) at p = that vectors ffc are linearly dependent 
over C ■ 

Proof of Theorem [SJ From Lemma 2 it follows that the rank of matrix F in flHI?|) is n, 
thus vector columns of F and G span the same linear space and therefore 

n 

Gk = Y,F^Ak (41) 

i=l 

or 

Dl{g,) = J2DmAk (42) 

Substitution of fj^Tl) in flH^ gives F(v4^^ + A)Ft^ = and since rankF = n we get 
A^^ + A = 0. From iMh it follows that 



i=l i=l 

which leads to Dx{Ai]^) = (since rankF = n). M 
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